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We investigate optically detected magnetic resonance signals from an ensemble of nitrogen- vacancy 
centers in diamond. The signals are measured for different light powers and microwave powers, and 
the contrast and linewidth of the magnetic-resonance signals are extracted. For a wide range of 
experimental settings of the microwave and light powers, the linewidth decreases with increasing 
light power, and more than a factor of two "light narrowing" is observed. Furthermore, we identify 
that spin-spin interaction between nitrogen-vacancy centers and substitutional nitrogen atoms in 
the diamond leads to changes in the lineshape and the linewidth of the optically detected magnetic- 
resonance signals. Finally, the importance of the light-narrowing effect for optimizing the sensitivity 
of magnetic field measurements is discussed. 

PACS numbers: 76.30.Mi, 76.70.Hb, 81.05.ug, 76.60.Es 



I. INTRODUCTION AND MOTIVATION 

A. Nitrogen-vacancy centers 

Nitrogen-vacancy (NV) centers in diamond have re- 
cently attracted attention within the fields of quan- 
tum information science^ ^ and optical magnetometrj/^ . 
Much of the focus has been on single NV centers, and 
impressive results have been obtained, such as quan- 
tum entanglement^ and nanoscale imaging with single 
spins^^. Here we focus on an ensemble of NV centers 
which can be utilized as a very sensitive /im to mm sized 
sensor for magnetic field&^E^. 



B. Optical magnet ometry 

In optical magnetometrjBSI^ magnetic resonance be- 
tween two ground states, which are split by a magnetic 
field is excited with a properly modulated light or 
magnetic field, where the modulation typically corre- 
sponds to radio- frequency or microwave (MW) range. 
The magnetic resonance is measured optically and the 
magnetic field is determined from the center modulation 
frequency corresponding to the magnetic resonance. 
To be able to determine the center frequency of the op- 
tically detected magnetic resonance (ODMR) signal and 
thereby the magnetic field with high precision, a narrow 
linewidth Au is desired. The magnetometer signal is pro- 
portional to C/Ai^, where C is the contrast of the ODMR 
signal. If the main source of noise is shot noise of the 
detected light, the magnetic field-sensitivity (measured 
in T/a/Hz) and also the noise-to-signal ratio will be pro- 
portional to Au/ ^Ca/P^, where P is the light power. 

Both the contrast and the linewidth depend on the MW 
power and the light power. To achieve high magnetic 
field-sensitivity one therefore needs to optimize the MW 
power and the light power. 



C. Light narrowing 

The dependence of the magnetometer signal on MW 
power and light power is now discussed. In general, both 
the contrast and the linewidth increase with MW power. 
In atom-light interactions, typically, the linewidth also 
increases with increasing light power. However, in some 
situations, the linewidth can decrease with increasing 
light power. This phenomenon is called light narrowing, 
with an early observation in an optically pumped alkali 
vapor reported in Ref. JJj In that work power broad- 
ening due to the MW field was reduced when the light 
power was increased. A similar light-narrowing effect on 
NV centers in diamond is the topic of this paper. A dif- 
ferent type of light narrowing, which occurs even in the 
absence of MW power broadening, was demonstrated in 
Ref. 12 and 13, There the linewidth was reduced due 
to the suppression of spin-exchange collisions when the 
alkali vapor became more and more polarized with in- 
creasing light power. This latter effect, however, does 
not appear to have a direct relation to this work. 



II. THEORY 

A. The Bloch formula 

Light narrowing of a power-broadened magnetic reso- 
nance can be understood using a two- level model pictured 
in Fig. [TJa). Consider a single NV center (or atom) with 
two relevant ground states |0) and |1) which are coupled 
by a MW magnetic field with coupling strength given by 
the on-resonance Rabi frequency Qr oc VPmw, where 
is measured in rad/s and Puw is the MW power. 
The longitudinal and transverse relaxation times of the 
atom are denoted Ti and T2 , and the corresponding decay 
rates are denoted 71 = 1/Ti and 72 = I/T2. We assume 
that the system can be optically pumped into the state 
|0) with a "pump rate" Fp. For a continuous- wave MW 
magnetic field with detuning A (in rad/s), one calcu- 
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lates using the Bloch equations (Appendix A) that the 
magnetic-resonance signal will be Lorentzian as a func- 
tion of MW detuning [Eq. (A2)]. After converting from 
rad/s to Hz, we find the signal as a function of frequency 
V in Hz to be 



S{u - uq) = S{oo) 



C • {Aiy/2f 



(1) 



5'(oo) is the signal when the MW's are off-resonant 
(|i^ — i^ol co), C is the contrast of the resonance, and 
Ai^ is the full width at half maximum (FWHM) measured 
in Hz. In our experiment, the signal is the detected flu- 
orescence from the NV centers. The signal is maximal 
off-resonance because optical pumping of the NV center 
leads to increased population of the |0) state which fluo- 
resces stronger than the |1) state. When the MW field is 
on-resonance, population is transfered from the |0) state 
to the |1) state leading to reduced fluorescence. 

The width of the resonance can be expressed as 



(2) 



where = VLr/ (27r) is the Rabi frequency measured 
in Hz. The "effective" relaxation times depend on the 
optical-pumping rate: l/Tf = 1/Ti + Ep and l/Tf = 
I/T2 + Fp/2. For low MW power the width of the mag- 
netic resonance is approximately Au ^ 1/ (ttTI^). In- 
creasing the MW power leads to power broadening, and 
for high MW power the FWHM 



Tf/Tf • /p. 



(3) 



equals the Rabi frequency magnified by the factor 2 • 
^Tf/Tf. 

For an ensemble of NV centers in type lb diamond 
(nitrogen density ~ 100 ppm), representative values for 
room-temperature intrinsic relaxation times (i.e., in the 
absence of light and microwaves) are TJ"*^ ^ 1 ms and 
T2"*^ ~ 1 /is and thereby differ by 3 orders of magni- 
tude. In our experiment, laser light is used to probe the 
NV centers and simultaneously optically pump the NV 
centers into the state |0). This optical pumping reduces 
Tf^. In general, laser light also leads to a reduced Tf^, 
but for certain "low" light powers, we can be in a regime 
where Tf^ is almost unaffected by the light but Tf^ is 
significantly reduced. For a power broadened magnetic- 
resonance line, the width will therefore decrease with in- 
creasing light power since the width is proportional to 



are situated in different local environments. This inho- 
mogeneous broadening can be modelled as though the 
NV |0) ^ |1) resonance frequency has a probability 
distribution P {1^0) [the probability per unit frequency 
for the resonance frequency to be at this value] with 
FWHM Ai/inh oc I/T2*. In general T2 < 2Ti and typically 
2^mtr ^ solids. The ensemble averaged ODMR sig- 

nal is then 



(4) 



where S{iy — uq) is the ODMR signal [Eq. ([T])] for a single 
NV center with resonance frequency uq. The FWHM 
Autot of the ensemble averaged ODMR signal can be 
found by calculating the FWHM of 5'tot(^) in Eq. Q. 
If the resonance frequency has a Gaussian probability 
distribution: P(i^o) oc exp (— z/Q/2cr^) [where a is defined 
such that the FWHM Az/inh ^ 2.35(j], the lineshape of 
Stoti^) will be the Voigt lineshap^^ which can be cal- 
culated using complex error functions. If on the other 
hand, the resonance frequency has a Lorentzian proba- 
bility distribution, the lineshape of 6'tot(^) will also be 
Lorentzian with a FWHM given by 



Ai^tot = Aui 



inh 



(5) 



where Aui^ih = 1/ (7rT2*) and Au is given by Eq. ([2]). Al- 
though it is probably most natural to assume that P(i^o) 
is Gaussian, we will assume that P{i^o) is Lorentzian. 
This is done mainly because one can obtain a simple an- 
alytical formula for the total width [Eq. ([5|], which later 
can be used for fitting the experimental data. Also, it 
seems that the measured ODMR spectra can be fitted 
well with Lorentzian functions (see for instance Fig.[2]3). 

Inserting the Bloch formula given by Eq. ^ into Eq. 
^ we find 




where we first assumed that the light does not affect the 
transverse relaxation rate 7!^ ~ 72"*^, and next omitted 
the term 72"*^ /tt because it is much smaller than Az/inh 
and will be hard to estimate from the width of the inho- 
mogeneously broadened ODMR lineshape. 



III. EXPERIMENT 



NV center level structure 



B. Inhomogeneous broadening 

Equation ^ is valid for a single two-level atom or 
for an ensemble of identical two-level atoms. However, 
the NV centers in the ensemble are different since they 



NV centers in diamond occur in at least two charge 
states, NV^ and NV~, where the negatively charged cen- 
ter is of primary interest since it is possible to optically 
polarize and readout the NV~ spin state. Figure [ijb) 
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FIG. 1: (a) Two- level model of ODMR dynamics. Two states 
|0) and |1) are coupled with a MW field (Rabi frequency Qr 
and detuning A); Fp optical-pumping rate; 71 longitudinal 
spin-relaxation rate; 72 transverse spin-relaxation rate, (b) 
Level scheme of the NV~ center, (c) Schematic of the confocal 
setup, (d) Points: Measured red fluorescence Pfi as a function 
of green pump power P. Solid line is a fit to a function of the 

ioim Pfl=k-P/{1 + P/Psat). 



shows the level structure of the NV~ center. The ground 
state is a spin-triplet and the states \ms = 0) and 
\ms = ±1) (the quantum number is the electron-spin 
projection along the NV axis) have a zero-magnetic-field 
splitting of 2.87 GHz. Applying a static magnetic field 
along the NV axis leads to a splitting of the \ms = +1) 
and \ms = —1) states. The ground states are coupled 
to a spin-triplet excited state using laser light. The 
zero-phonon line (ZPL) occurs at 637 nm, however, typ- 
ically the NV center is excited with green light (532 nm) 
through the phonon sideband (FSB). Once excited, an 
NV~ center can decay back to the ground states directly 
(through the spin-conserving transition ^ ^E) emit- 
ting red fluorescence in the FSB (with wavelength ^ 
637-800 nm) or through a spin non-conserv ing c hannel 
through the spin-singlet states ^Ai and ^E l^'^l^^l (The 
order of the singlet levels shown inFigure[TJb) is given in 
accordance with the recent woiW^^.) The latter chan- 
nel through the singlet levels is more probable for the 
\ms = ±1) states than for the \ms = 0) state of^E. This 
leads to the possibility to efficiently optically pump, po- 
larize and read-out the NV~ spin state. 



B. Setup and procedure 

The setup is shown in Fig.[TJc). We use the high pres- 
sure high temperature HFHT diamond sample S5 (inves- 
tigated previously in Ref. |7j) with the following proper- 
ties: nitrogen concentration ^ 200 ppm, NV~ concen- 
tration ^ 12 ppm, and a cut perpendicular to the [100] 
crystallographic direction. The diamond is located in 
a static magnetic field produced by a permanent mag- 



net. A MW generator sending current through a wire 
located on top of the diamond produces the MW field. 
Green (532 nm) light from a laser is focused onto the di- 
amond using a confocal setup with an aspheric lens with 
focal length of 4 mm and numerical aperture NA=0.60, 
and the red flouresence is detected with a photodiode 
after passing through a dichroic mirror. The amount 
of red fluorescence as a function of green light power 
is plotted in Fig. [TJd) together with a fit of the form 
Pfi = k- P/(l + P/Psat) where k = 6.21(2) • 10"^ and 
Psat = 4.8(2) W are fitted parameters. The saturation 
power for equalizing the populations of the ground (^^2) 
and excited {^E) states is denoted Psat- For the light 
powers investigated here (P < 500 mW), we hardly see 
any saturation. The MW frequency is scanned and the 
ODMR signal is recorded using an oscilloscope. Four 
different orientations of the NV axis in the diamond lat- 
tice are possible leading to eight magnetic resonances 
(since each of the four orientations has two electron- 
spin resonances: \ms = 0) ^ \ms = +1) and \ms = 0) ^ 
|ms = — 1)). A static magnetic field of 77 G is aligned 
along one of the [111] crystallographic directions. In this 
case, NV centers oriented along this axis have magnetic- 
resonance frequencies which are different from the over- 
lapping resonances of the three other orientations. For 
this magnetic-field direction, one should therefore ob- 
serve four magnetic-resonance frequencies. 



C. Optically detected magnetic resonance signals 

ODMR signals were recorded for different settings of 
the MW and light powers. To determine the strength of 
the MW field, we employed a pulsed calibration measure- 
ment where the on-resonance Rabi frequency was mea- 
sured for different MW settings. We found the expected 
fn oc ^JPmw square-root dependence of the Rabi fre- 
quency on the MW power. 

Examples of ODMR signals are shown in Fig. |2] The 
ODMR signals are normalized to the signal obtained 
when the MW frequency is far detuned. Figure [2ja) 
shows an ODMR signal for a wide scan range of the MW 
frequency. Four magnetic-resonance frequencies of the 
NV center are observed, as expected when the magnetic 
field is along one of the [111] directions. We now focus 
on the magnetic resonance with the lowest resonance fre- 
quency, which corresponds to the Im^ = 0) ^ Im^ = — 1) 
transition of NV centers oriented along the magnetic field 
direction. A narrower scan of the lowest resonance is 
shown in Fig. ^^). The MW power was set at a low 
value for this measurement (corresponding to Jr = 0.02 
MHz). In this case, the resonance is sufficiently nar- 
row to partially resolve the hyperfine structure of the 
NV center. [The hyperfine interaction of the NV elec- 
tron spin and the ^^A^ nuclear spin 1 = 1 splits each 
\ms = 0) ^ |ms = ±1) resonance into three resonances 
with splitting Ahf = 2.2 MHz (see for instance Fig. 10 in 
Ref. 7).] 
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FIG. 2: Examples of ODMR signals taken with the S5 (a)-(e) and S3 (f) diamond samples at rooom temperature, (a) ODMR 
signal with a wide scan range of the MW frequency, (b) Narrow scan range ODMR signal of the magnetic resonance centered 
around 2654 MHz. (c) Light power fixed and four different settings of the MW power, (d) MW power fixed at a relatively 
high value and four different settings of the light power, (e) Demonstration of light narrowing with the S5 sample. Notice the 
two y-axes. (f) Demonstration of light narrowing with the S3 sample. The bottom spectrum (P = 500 mW setting) has been 
displaced vertically by -0.001. 



Figure 2jc) shows ODMR signals for different settings 
of the MW power. The central resonance frequency 
uq ^ 2654 MHz shifts slightly to lower frequencies with 
increasing MW power due to heating^^. Also, one clearly 
sees that increasing the MW power leads to power broad- 
ening of the ODMR signal. 

We now proceed with the demonstration of light nar- 
rowing. Figure [2|d) shows examples of ODMR signals for 
four different settings of the light power P ranging from 
0.02 to 500 mW. The MW Rabi frequency was at a rela- 
tively high value of 1.1 MHz. Fitting the spectra reveals 
that the linewidth decreases with increasing light power 
(the details of the fit procedure and the fit results are 
discussed later). This is the light-narrowing effect. The 
spectra with the lowest and highest light-power settings 
are also plotted in Fig. [2|e) [notice the two y-axes in Fig. 
[2|e)]. There we clearly see a reduction of the linewidth 
by approximately a factor of two when the light power is 
increased from 0.02 mW to 500 mW. 

We also demonstrate light narrowing with the S3 di- 
amond sample which has a lower concentration of both 
NV~ centers (~ 0.012 ppm) and nitrogen impurities (< 1 
ppm) (for details of the S3 sample see Ref. 7). Due to 
the lower impurity concentration, the S3 sample has less 
inhomogeneous broadening compared to the S5 sample 
(T2* ^ 300 ns for the S3 sample compared to ^ 100 
ns for the S5 sample). This makes the hyperfine struc- 
ture better resolved in the S3 sample compared to the 



S5 sample. Figure [2jf) shows ODMR spectra for the S3 
sample for two settings of the light power. Notice that 
the bottom spectrum (P = 500 mW setting) has been 
displaced vertically by -0.001 for better visualization. At 
the low light power setting, the spectrum is relatively 
broad and a single resonance is observed. At the high 
light power setting, the spectrum is narrower, and the 
hyperfine structure can be partially resolved. Le., spec- 
tral features which are not resolved at low light power 
can be resolved at higher light power due to the light- 
narrowing effect. 



D. Side resonances 

A pair of side resonances at a distance of ~ 33 
MHz from the central resonance can be seen in Fig. [2|c) 
(marked by arrows). The contrasts of the side resonances 
increase with increasing Rabi frequency. The side reso- 
nances are due to MW-induced simultaneous spin-flips 
of NV centers and nearby substitutional nitrogen atoms 
(also called PI centers)P'. 

Let be a MW spin-transition frequency of the NV 
center and S a transition frequency between two spin 
states of a nearby PI center (which has electron spin 
1/2 and nuclear spin 1). Due to magnetic dipole-dipole 
interaction between the NV center and the nearby PI 
center, the spin eigenstates of the NV center and the PI 
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center will be slightly mixed. In this case, it is possible to 
drive MW transitions with resonance frequencies z^o i ^ 
which flip the spin of the NV center and the PI center 
simultaneously. 

The side resonances were studied in detail and the 
results will be reported elsewhere^^. In short, many 
side resonances were observed, with resonance frequen- 
cies depending on the magnetic-field amplitude and 
orientatioiP^. The contrasts of the side resonances were 
found to be sample dependent and to increase linearly 
with MW power at low MW powers and saturate at 
higher MW power4^. In the S5 sample [see Fig. [2|c)], 
the side resonances were rather small and did not over- 
lap much with the resonance centered at , but in other 
samples the side resonances were larger and in some 
cases they significantly distorted the spectrum of the NV 
centei^^. We note that while side resonances were ob- 
served in the S5 sample, no side resonances were observed 
in the S3 sample. This makes sense since the concentra- 
tion of PI centers in the S3 sample (< 1 ppm) is much 
lower than in the S5 sample 200 ppm)P. We would 
therefore not expect any visible side resonances in the S3 
sample. 

The rate of MW-induced simultaneous spin-flips 
7NV-P1 should depend on several parameters including 
PI concentration, MW detuning from the side resonances 
at uq ± MW power and light power. We note that the 
mechanism behind the MW-induced simultaneous spin- 
flips is rather complicated. A detailed analysis of how 
NV-Pl interactions change the ODMR spectrum is be- 
yond the scope of this paper. We would instead like 
a simple model for how the MW-induced simultaneous 
spin- flips of NV and PI centers influence the width of 
the magnetic resonance centered at uq. Spin- flips of the 
NV electron spin due to interaction with PI centers can 
be thought of as a source of Ti -relaxation for the NV cen- 
ters. We will assume that the only effect of NV-Pl inter- 
action is to increase the Ti -relaxation rate by an amount 
corresponding to the spin-flip rate: 71 = 7}"*^ +7nv-pi- 

For simplicity, we will neglect the detuning dependence 
of the spin-flip rate. This assumption should be valid 
since the side resonances are well separated from the res- 
onance centered at and since we only are interested 
in the effect of the simultaneous spin-flips on the width 
of the resonance centered at uq. 

Since the contrast of the side resonances increase with 
increasing MW power, the spin-flip rate should also in- 
crease with increasing MW power. The exact dependence 
on the MW power is not known, but we will assume that 
the rate is proportional to MW power, or equivalent, the 
square of the Rabi frequency. We will also assume that 
the rate of spin-flips saturate for high MW powers. We 
end up with the following empirical formula for the spin- 
flip rate: 

7NV-P1 = a{P) ■ fll (1 + fl/Q (7) 

The saturation Rabi frequency /o and the proportionality 
constant a{P) and its dependence on light power will be 



determined from the experimental data. 



E. Fitting procedure 

The normalized ODMR signals S{v) are fitted to the 
following function 



mj=-i - ^0 - m/ • Ahf) +r 

which is a sum of three Lorentzians each with peak height 
(or amplitude) A, center frequency + mj • Ahf (m/ = 
-1,0,1) and FWHM 2g. This fit function with three 
components separated by Ahf = 2.2 MHz is chosen due 
to the hyperfine structure of the NV center. 

The fits are shown on top of the data in Fig. [2] as 
solid lines. From the fits we extract the FWHM 2g and 
the amplitude A of the ODMR signals. Notice that we 
choose to define the width of the ODMR signal as the 
width of an individual hyperfine component (which is 
smaller than the width of the total ODMR signal). When 
fitting, we excluded the part of the spectrum where the 
side resonances are located (~ 20 MHz of MW scan range 
centered at each side resonance was excluded). Fitting 
the whole ODMR spectrum did not change the fitted 
parameters significantly. 



F. Linewidth 

More ODMR signals were recorded, and the fitted 
FWHM and amplitude are plotted in Fig. |3] as a func- 
tion of light power and Rabi frequency, respectively. The 
light power was varied more than four orders of magni- 
tude and we therefore use logaritmic scale for the power 
axis. 

We will first focus on the dependence of the linewidth 
on the light power. Figure |3jb) shows the width as a 
function of the light power for several settings of the 
Rabi frequency. We observe a significant light narrow- 
ing (decreasing width with increasing light power) for all 
settings of the Rabi frequency. A maximum of 2.5 reduc- 
tion of the width was found for the higher Rabi frequen- 
cies. Note that in a related stud}'^^, where the widths of 
ODMR signals from single NV's were examined, no light 
narrowing was observed. In that stud}P^ the light inten- 
sity was high and lead to significant power broadening 
instead of light narrowing. 

For the highest light power (P = 500 mW), the width 
seems to have a linear dependence on the Rabi frequency 
[Fig. [sja)]. For lower light powers, we observe a non- 
linear scaling with Rabi frequency, in particular we found 
that the width as a function of Rabi frequency has a neg- 
ative curvature (i.e., the second derivative of the width 
with respect to Rabi frequency is negative) . According to 
Eq. ([6| for the linewidth, the curvature is negative if 71 
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FIG. 3: Points: Experimental results of the FWHM and amplitude of ODMR signals as a function of Rabi frequency and light 
power. Solid lines: Global fits to the FWHM and the amplitude. The uncertainties on the experimental data points correspond 
to the standard deviations found from repeated measurements of the width and contrast. 



increases with Rabi frequency. We believe that the neg- 
ative curvature is a consequence of MW-induced simul- 
taneous spin- flips of NV and PI centers which effectively 
increases 71 with increasing Rabi frequency. 

As a remark, we note that the presence of side reso- 
nances was correlated with a negative curvature. In the 
S5 sample (which has a high density of PI centers), both 
side resonances and a negative curvature were observed, 
while in the S3 sample (which has a low density of PI 
centers), neither side resonances nor negative curvature 
were observed. This supports our hypothesis, that the 
negative curvature is due to NV-Pl simultaneous spin- 
flips. 

Assuming that the longitudinal spin-relaxation rate is 
7i = Tr*"" + 7NV-P1, where 7nv-pi is given bvEq. 
the linewidth of the ODMR resonance [Eq. ([6|] is ex- 
pected to be on the following form: 



Ai/tot ~ Az^inh + fn • 



4^mtr 



(1 + W) 



7r^^ + a(P)./H/(l + /I//o^) + c-p- 



(9) 

For low light powers, the optical pumping rate is assumed 
to be proportional to light power: Tp = cP with propor- 
tionality constant c. The pump rate should saturate at 
higher light powers due to the finite lifetime of the singlet 



states responsible for the optical-pumping mechanism of 
the NV center. The effect of saturation on the width of 
the ODMR spectrum is analyzed in Appendix B using a 
five-level model for ODMR dynamics. We find that one 
should include the term (I + P/Pq) in Eq. (|9| for the 
width [see Eq. (B4) in Appendix B]. Note that the opti- 
cal pumping saturation power Pq which is related to the 
singlet state lifetime (^ 200 ns at room temperatur^^ 
should be smaller than the saturation power Psat [shown 
in Fig. [ijd)] for equalizing the populations in the ground 
and excited states {Psat is related to the excited state 
lifetime ^12 ns). 

We choose to make a global fit of the width as a func- 
tion of Rabi frequency and light power to Eq. (|9|) with the 



fit parameters: { Ai/inh, 7^*77^"*', a/Y^''^', cTVf *^ Pq}. 
Since data with 12 different light-power settings were fit- 
ted, the fit parameter a(P) 772"*^ will be a vector with 
12 components. The fits to the experimentally measured 
widths are shown in Fig. [sj^a) with solid lines and we see 
that the measured data can be well fitted by Eq. The 
values of the fit parameters are given in Table |I] and Fig. 
Q. We note that the uncertainty on the fit parameter 
a(P)/Y2^^^ is large for high light powers because the width 
depends only weakly on 7nv-pi when cP ^ 7nv-pi- 
We emphasize that the term a- f'^/ (l + /^//o ) in Eq. 
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(|9| is necessary to obtain good agreement between data 
and fit since this term is responsible for the negative cur- 
vature of the width as a function of Rabi frequency. The 
fitted parameter is shown in Fig. Q and is found 

to be in the range (0.1 — 0.3)/MHz for the different light- 
power settings. 

The term a • /d/ (l + fn/fo) the fit function for 
the width [Eq. ^] is interpreted as the relaxation rate 
7NV-P1 due to NV-Pl interaction. For a Rabi frequency 
of 1 MHz, fit parameters ^ 0.15/MHz and /o = 

1.0 MHz and an estimated 72"*^ = l//is, we calculate 

7NV-P1 = Ih""'' ■ {a/lf'') ■ fn (1 + /I//0) « 0.0757r'' « 
0.075//is. For a Rabi frequency of 1 MHz, the contri- 
bution from 7NV-P1 to the effective Ti-relaxation rate 
7®^ is therefore much larger than the contribution from 
the intrinsic relaxation rate which is estimated to be 
^mtr ^ 1/ (1 ms) = 0.001//is. The intrinsic longitudi- 
nal spin-relaxation rate is therefore negligible compared 
to the NV-Pl spin flip rate at high MW powers. We 
conclude that intrinsic Ti -relaxation will probably not 
be important for continuous- wave (CW) magnetometry 
using NV centers in diamond samples with comparable 
nitrogen concentration. 

We choose to fit the obtained a(P)/72"*^ to an empiri- 
cal function 



aiP/{l + P/bi)+ci. 



(10) 



The fit parameters are given in Table |T] and the fit funtion 
plotted in Fig. Q as a solid line fits the data well. It is 
not clear why the parameter a should have this particular 
dependence on light power. In order to better understand 
the process of MW-induced NV-Pl simultaneous spin- 
flips and the associated spin-flip rate, one would need 
to further develop a theoretical model and probably do 
more direct measurements of the NV-Pl spin- flip rate 
instead of these indirect measurements where the NV-Pl 
spin flip rate is inferred from the width of the ODMR 
spectra. 

We also note that 71^*77^"*' was fitted to 0.0014(3). 
This agrees well with the expected value ^ 10~^, which 
was calculated using estimated intrinsic relaxation rates 
(^^ntr ^ i/ms and 72^*^ ~ l/jus). We can also estimate 
the value of cPq from the fit parameters. We have cPq = 
72''*'- -^0 ~ 7^"*'-0.018/mW-39mW ^ 0.77^"*^ ^ 

0.7//is. As calculated in Appendix B using a five-level 
model for ODMR dynamics, cPq should equal Fg, where 
Fg ^ 1/ (200ns) = 5//is is the decay rate from the singlet 
state. The value of cPq estimated from the fit parameters 
is therefore a bit on the low side compared to the value 
estimated from the singlet state lifetime. 



G. Contrast 

In order to have high magnetic sensitivity, the ODMR 
resonance should both be narrow and have a high con- 
trast. The contrast C is defined as the total peak height 
of the ODMR signal: C = 1 - S{iyo)/S{oo) [see Eq. Q], 
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FIG. 4: Points: Fit parameter a/Y^^^ for different light pow- 
ers. The uncertainties on the points equals the 68% confi- 
dence intervals for the fit parameters a[P)/Y^^^ calculated 
from the non-linear least squared fit of the measured widths 
to Eq. Line: Weighted fit of a(P)/7^^*' to the function 
aiP/(l + P/6i)+ci . 



Az^inh 


3.08(9) MHz 


^intr/^intr 


0.0014(3) 




see Fig. (Ul 




0.018(5) 1/mW 


Po 


39(13) mW 


fo 


1.0(1) MHz 


ai 


0.5(2) l/(MHz • mW) 


bi 


0.5(2) mW 


Cl 


0.074(7) 1/MHz 





22.9(4) • 10"^ 




0.71(6) mW 




0.0047(6) (/is)"^ 



TABLE L Fit parameters. 1st column: fit parameters. 2nd 
column: values for the fit parameters. The fit functions are 
given by Eqs. ([9|, ([lO| and (11). 



which is larger than the amplitude or peak height A of 
an individual hyperfine resonance [see Eq. (|8|]. We can 
calculate the relation C = A • [l + 2^^/ [A\^ + g'^)] from 
Eqs. ([1]) and (|8|. When the width is much larger than the 
hyperfine splitting {g^ ^ A'^^) we find a simple relation 
C = 3A. 

Figure [sj^d) shows the fitted amplitudes as a function 
of light power for several different MW settings. The 
measured amplitudes lie in the range (0.1 — 6.8) • 10~^. 
The maximum measured contrast is therefore C ^ 3 • 6.8 • 
IQ-^ ^ 2.0%. As seen in both Fig. ^c) and Fig. [sjc), 
the contrast increases with MW power. The dependence 
on light power is more complicated [Fig. |3|d)]. First the 
amplitude increases with light power, then a maximum 
is reached, and finally the amplitude decreases. Both 
the maximum amplitude and the light power needed to 
achieve the maximum amplitude depend on the MW 
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power. 

The theoretical expected contrast was calculated using 
the two-level model described in detail in Appendix A. 
Upon rearranging, we obtain 



cP 



cP + 7™*'- (1 - e) 



(2^)- 



cP 



(11) 



Here is an overall proportionality constant defined in 
Appendix A. The factor 1/4 is included to take into ac- 
count that only 1/4 of the NV centers (the ones which 
are oriented in the direction of the magnetic field) are 
in resonance with the MW field. In Eq. ( pT] ) for the 
contrast, the pump rate is assumed to be Fp = cP. 
From Eq. (11), we see that the pump rate Fp should 



be larger than the Ti -relaxation rate 7^"*^ in order to 
achieve high contrast [Fp > 7^"*^ (1 - 6>)]. If the light 
power is too low, it will not be possible to polarize the 
NV centers. This explains the initial increase in contrast 
with increasing li^ht power. However, at high light pow- 
ers [Fp > (27r/p) /72"*^], the optical pumping becomes 
too fast compared to the MW driving field, and it will 
not be possible to create a population difference between 
the two ground states |0) and |1) involved in the MW 
transition. This explains the decrease in contrast at high 
light powers. 

We choose to make a global fit of the contrast as a 
function of Rabi frequency and light power to Eq. (11) 
with the fit parameters 
to the data are shown in 



,7i 



7^^*^}. The fit 



?,7l-*7c. 

Fig. [3] (c), (d) as solid lines 
and the fit parameters are given in Table |lj The fit func- 
tion Eq. (11) agrees reasonably well with the data. The 



discrepancy between data and fit is probably due to the 
simplicity of the two-level model leading to Eq. (11) for 
the contrast, where we did not include the effect of inho- 
mogeneous broadening, MW-induced simultaneous spin- 
fiips of NV and PI centers, and possible population of 
the singlet and excited states. 



H. Magnetometer sensitivity 

Using a measured ODMR spectrum, one can calculate 
the experimental magnetometer sensitivity from the un- 
certainty on the fitted center frequency and the averaging 
time used to obtain the ODMR spectrum. However, since 
our measurement uncertainty is mainly limited by tech- 
nical noise specific to our apparatus, this will not tell us 
about the fundamental limits to the sensitivity. We in- 
stead investigate the projected sensitivity for a shot-noise 
limited measurement. 

The shot-noise limited sensitivity is of the order of 



Sb 



27r Az/ 



where Az/ is the FWHM of the magnetic resonance, 
7 = 1.761 • 10^^5~^T~^ is the gyromagnetic ratio for the 
NV center and IZ is the rate of detected photons. The 



(12) 



experimental sensitivity Sb is calculated using Eq. (12) 
and experimentally measured values of Az^, C and IZ. 
When calculating Sb-, Az^ was set to the fitted FWHM, 
C was set to 3A (these assumptions are valid when the 
hyperfine structure is not resolved) and IZ was estimated 
from the measured power of the fiuorescence. We can 
calculate IZ = Pji/Efi where Eji = hc/X {h Planck's 
constant, c speed of light, A ~ 670 nm wavelength of the 
fiuorescence). Using the calibration shown in Fig. [TJd), 
we find that 1 mW of green light produces 6.2 /iW of 
red fiuorescence which corresponds to a rate of detected 
photons 7^ = 2.1 • 10^^ s-^ 

Figure |5] shows a contour plot of 5'p as a function of 
P and /p. In total, experimental data with more than 
100 different settings of (P, /p) were used to create the 
contour plot. The best sensitivity Sb ^ 0.1 nT/V^Hz is 
reached for the highest light power P = 500 mW and 
Rabi frequency settings /p = 0.34, 0.62 and 1.14 MHz. 
For even higher light powers (P > 500 mW), the sensitiv- 
ity is expected to eventually become wors^^. We note 
that our best projected shot-noise limited sensitivity is 
significantly better than the actual measured sensitivity 
obtained with ensembles of NV centers so fai^ 

When detecting fiuorescence from a diamond with a 
confocal microscope, the photon detection efficiency is 
rather small (in our case ~ 0.8 %) due to the high index 
of refraction of diamond (n ~ 2.4). We note that there 
exist several methods to increase the detection efficiency 
including the "side-collection" technique^, the use of solid 
immersion lenses^^ and a detection scheme based on 
IR absorption^. If such techniques would be employed in 
our experiment, the detection efficiency and thereby the 
rate of detected photons could be increased by roughly 
two orders of magnitude. This would lead to an over- 
all improvement of our best projected shot-noise limited 
sensitivity by a factor of 10 to the value 10 pT/>/Hz. 

Combining the fact that the settings of (P, /p) where 
the sensitivity is optimal [Fig. [5], are also the settings 
where the width of the magnetic resonance is reduced by 
the light-narrowing effect [Fig.[3];b)], we reach the con- 
clusion that the light-narrowing effect plays an important 
role for optimal sensitive magnetometry with NV centers 
in diamond. 



I. Light narrowing and IR absorption 

We have so far focused on ODMR signals detected us- 
ing ffuorescence from the spin-triplet excited states. An 
alternative way of measuring magnetic fields using NV 
centers in diamond is to use infrar ed (I R) absorption 
at the singlet transition at 1042 nnF*^. In this case, 
a green laser is used to pump NV centers into the lower 
spin-singlet state, the MW frequency is scanned around 
magnetic resonance while the absorption of an IR laser 
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sensitivity [nT/VHz] 




10 10 10 10 

light power [mW] 

FIG. 5: Contour plot of the experimental magnetic field- 
sensitivity Sb SiS a function of light power P and Rabi fre- 
quency fn. Notice the logarithmic scales. 



(1042 nm) is measured. 

We have analyzed a five-level model for the NV cen- 
ter. The details of the model are given in Appendix 
B. As shown in Fig. |6| the model includes two ground 
states |0) , |1), two excited states |eO) , |el) and one state 
\s) representing the singlet states. The IR absorption is 
proportional to the number of NV centers in the singlet 
state aiR oc pss (here pss denotes the diagonal density- 
matrix element associated with the state \s)). Using 
Bloch equations (see Appendix B), one can calculate that 
the steady-state solution pss(A), where A is the MW de- 
tuning, will have the same lineshape and the same width 
as the ODMR signal S{A) measured using fluorescence. 
Light narrowing should therefore also occur when mea- 
suring IR absorption (in the sense that the width of the 
ODMR signal measured using IR absorption decreases 
with increasing green pump power). This makes sense 
since light narrowing of the MW resonance should not 
depend on the detection method (fluoresence vs absorp- 
tion). 



IV. CONCLUSIONS 

In this paper, the light-narrowing effect was discussed 
using a two-level model, which is sufficient to demon- 
strate the basic phenomenon. We experimentally demon- 
strated light narrowing using NV centers in diamond and 
found a reduction in the linewidth of the microwave tran- 
sition by more than a factor of two when the light power 
was increased. The obtained reduction in linewidth is 
compara ble to what was found in earlier studies of light 
narrowin^^^ in alkali vapor cells. It was also found that 
MW-induced simultaneous spin-flips of NV centers and 
PI centers are important for determining the lineshape 
and width of the ODMR signal. Finally, we showed that 
light narrowing is an essential factor in optimizing the 



sensitivity of magnetometers based on NV centers in di- 
amond. 
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Appendix A: Two-level model 

We analyze the situation of a two-level atom driven 
by a MW field [see Fig.jlja)] with detuning A and power 
corresponding to the on-resonance Rabi frequency Qr 
A and are defined to be angular frequencies measured 
in rad/s. We assume that the upper level can be opti- 
cally pumped into the lower level with a rate Fp. The 
populations of the |0) and |1) states are assumed to de- 
cay towards equal populations at a rate 71 = 1/Ti. The 
transverse relaxation rate is denoted 72 = I/T2. We have 
the following equations of motion (the Bloch equations) 
for the density-matrix elements: 



Poo 


= - — (Poi - 


71 

p\o) - — {poo - 


pll) + Fppii 


(Ala) 


pll 


= +^(poi- 


71 

Pio) - y {pll - 


Poo) - Fppii 


(Alb) 


Poi 




^ .^R . 

Poi (Pii 


- Poo) 


(Ale) 


pio 


= -(7f + ^A) 


.^R . 
pio - [pll 


- poo) . 


(Aid) 



The effective dephasing rate is denoted 7!^ = 72 + Fp/2. 
The s t eady-s tate solutions can be calculated from Eqs. 
(Ala)-(Ald). Similar to what was done in Ref. 24^ we 
assume that the ODMR signal can be written as 5'(A) = 
<^Poo + /^Pii7 where a > /3. In this case, the ODMR signal 
will have the following form: 



5(A) = S{oo) 



1 



C7^ 



A2+72 



(A2) 



S{oo) is the signal when the MW's are off-resonant 
(|A| 00) and 7 is the hwhm in rad/s. We calculate 
the FWHM in Hz: Ai/ = 7/7r to be 



Az/ = 




47f 



(A3) 



where fn = 0,^,/ (27r) is the Rabi frequency in Hz. The 
ODMR signal contrast is 



c = e- 



mi 



(A4) 



where = {a — /3) / (2a). Note that similar expressions 
for the linewidth and contrast have also been derived in 
the Appendix of Ref. |24i 
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FIVE LEVEL MODEL 



|el> 



|eO) 



|1> 
|0> 



A..; 



\s) 



y / ^ s 
/ 

/ 

/ 

■71 = 1/7^1 

72 = I/T2 



FIG. 6: Five-level model of ODMR dynamics. The two 
ground states |0) and |1) are coupled with a MW field (Rabi 
frequency and detuning A). Laser light excite the states 
|0) and |1) with equal rate Fp to the excited states |eO) and 
I el), respectively. The excited states decay to the ground 
state with rate Fq. The state |el) can also decay to the sin- 
glet state \s) with rate F/. The singlet state \s) decays with 
a rate F^ to either |0) or 



Appendix B: Five-level model 



We now present a five-level model describing light nar- 
rowing and optical pumping of the NV center. The model 
includes the ground states, excited states and singlet 
states of the NV center (see Fig. [6|, and is therefore 
more detailed than the two-level model [see Fig. [ija)] 
described earlier. Previously, a five level model has been 
used for calculating the dynamics of a sing le NV centeM 
We will use the five-level model to calculate the lineshape 
and width of ODMR spectra based on both fluoresence 
and infrared absorption detection. 

As in the two- level model, two ground states |0) and 
|1) are coupled by a MW field with Rabi frequency VLr 
and detuning A. Laser light can excite the ground states 
with a rate Fp through a spin-conserving transition to 
the excited states which in the model are denoted |eO) 
and I el). The excited states can decay with rates Fq 
to the ground states. In the model, the |el) state can 
also decay with rate F/ to the state \s) representing the 
singlet states. Finally, the state \s) can decay with a rate 
V s to the ground states |0) and |1) with equal probability. 
Realistic values of the decay rates are Fq = 1/(12 ns), 
F/ = 1/(12 ns) and F^ = 1/(200 ns). 

The equations of motion for the density-matrix ele- 
ments for the five-level system {|0) , |1) , |eO) , |el) , l^)} 



can be written as: 



Poo 



Poi 
Pio 

PeOeO 

Pelel 
9ss 



(poi - Pio) 



71 

2 



(p. 



'00 



OeO " 



+ ^ (Poi - Pio) - y (Pii - 
-(7f -iA)poi+i^(pii 



{^f + iA) pio 



Tppoo — ToPeOeO 

^PPW — ^OPelel — ^fPelel 
T/Pelel ~ ^sPss- 



Pll) 

(Bla) 

Poo) 

(Bib) 

-Poo) (Blc) 

-Poo) (Bid) 
(Ble) 
(Blf) 
(Big) 



The Ti -decay is here assumed to be towards equal pop- 
ulation in |0) and |1). The transverse relaxation rate 
is denoted 72 = I/T2, and the effective transverse relax- 
ation rate is defined as 7!^ = 72 + Fp/2. The NV centers 
are pumped on the spin-conserving transitions |0) |eO) 
and |1) I el) with off-resonant green 532 nm light. Co- 
herences between ground and excited states are set to 
zero since we use off-resonant excitation which does not 
preserve such coherences. 

The steady-state solutions for the above density- 
matrix elements can be calculated analytically from Eqs. 
(Bla)-(Blg). If one measures the amount of fluorescence 



from the spin-triplet excited states, the ODMR signal will 
be on the form 



^(A) (X PeOeO 



To 



/ 



-Pelel S{oo) 



A2 + 72 



(B2) 

where C is the contrast, 7 the hwhm in rad/s, and 6'(oo) 
is the obtained signal when the M W's are far off reso- 
nance. The lineshape given by Eq. ( |B2| ) is t he sa me as 
was obtained with the two-level model [Eq. ( |A2[ )]. We 
calculate the FWHM in Hz which under the assumptions 
Fq ^ Ej, Fp <C Fq and 71 <C F5 has the simple form 



if 



,eff 



i + rp/(4r,) 



71 + rp/4 



R- 



(B3) 



The assumption Fp <C Fq means that the light power 
should be sufficiently low such that the excited states are 
not significantly populated. This assumption should be 
valid in our experiments, since we do not observe much 
saturation of the amount of fluorescence with increasing 
light power [see Fig. [TJd)]. 

For low excitation power: Fp/4 <C F5, the two-level 
and the five- leve l models give identical results for the 
width [Eqs. ( |A3| ) and (B3)], if one associates the five- 
level excitation rate diveded by four with the two-level 
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optical pumping rate: f p/4 ~ Pp. For higher excitation 
powers: Vp/A > F^, the singlet states become partially 
occupied, leading to saturation effects. Thi s is included 
in the five- level formula for the width [Eq. (B3)] by the 



term 



i + rp/(4r,) 



Equation (|B3| can be re- written in terms of light power 



as 




47f -[l + P/Po] 
li+cP 



■fl (B4) 



In the above we use that the excitation rate is propor- 
tional to light power Fp = 4cP. The proportionality 
constant is denoted 4c. The optical pumping saturation 
power Po due to population of the singlet states is defined 
by the equation cPq = Tg. 



Instead of detecting fiuorescence emitted from the 
spin-triplet excited states, one can measure ODMR sig- 
nals using absorption at the infrared ( IR) ^ transition at 
1042 nm involving the spin-singlet states^ In this case, 
a green laser is used to pump NV centers into the lower 
spin-singlet state, the MW frequency is scanned around 
magnetic resonance while the absorption of an IR laser 
(1042 nm) is measured. The IR absorption will be pro- 
portional to the number of NV centers in the singlet state 
<^IR pss' One can calculate from Eqs. ( |Bla| )-(Blg) that 
the steady-state solution ^^^(A) will have th e sa me line- 
shape [Eq. ( |B2[ )] and the same width [Eq. (B3)] as the 
ODMR signal S{/S) measured using fiuorescence. Light 
narrowing can therefore also occur when measuring IR 
absorption (in the sense that the width of the ODMR 
signal measured using IR absorption decrease with in- 
creasing green pump power). 
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